Vortex lattice disorder in YBa 2 Cu 3 07_5 probed using /3-NMR 
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/3-detected NMR (/3-NMR) has been used to study vortex lattice disorder near the surface of the 
high- Tc superconductor YBa2Cu307-^ (YBCO). The magnetic field distribution from the vortex 
lattice was detected by implanting a low energy beam of highly polarized 8 Li + into a thin overlayer 
of silver on optimally doped, twinned and detwinned YBCO samples. The resonance in Ag broadens 
significantly below the transition temperature Tc as expected from the emerging field lines of the 
vortex lattice in YBCO. However, the lineshape is more symmetric and the dependence on the 
applied magnetic field is much weaker than expected from an ideal vortex lattice, indicating that 
the vortex density varies across the face of the sample, likely due to pinning at twin boundaries. At 
low temperatures the broadening from such disorder does not scale with the superfluid density. 

PACS numbers: 74.72.-h, 74.25. Qt, 75.60.-K, 75.70.Cn 



I. INTRODUCTION 

The vortex state of cuprate superconductors is of cen- 
tral importance in understanding high- Tc superconduc- 
tivity (HTSC). One the most well studied quantities is 
the internal magnetic field distribution p(B) associated 
with the vortex lattice (VL)ji~— As discussed below, sev- 
eral methods can be used to measure p{B), which de- 
pends on the London penetration depth A, the coherence 
length £ and, to a lesser extent, the internal structure 
of the vortices^ and non-linear and non-local effectsi^^ 
The form of p{B) has a distinctive asymmetric shape due 
to the spatial magnetic inhomogeneity characteristic of 
an ordered two-dimensional (2D) lattice of vortices. One 
basic feature in p(B) is a prominent high field tail asso- 
ciated with the vortex cores, which depends on the mag- 
nitude of £. There is also a saddle point in the local field 
profile located between two vortices. This gives rise to a 
Van Hove singularity or sharp peak in p{B) below the av- 
erage field. The overall width or second moment of p(B) 
depends primarily on the London penetration depth A, 
the lengthscale over which the magnetic field is screened. 
Anisotropy of the Fermi surface or the superconducting 
order parameter can result in a different VL but the main 
features are similar for any ordered latticei 10 ' 11 

Another general feature associated with any real VL 
is disorder arising from vortex pinning at structural de- 
fects and impurities . 12 ' 13 Structural defects are present 
in all superconductors to some degree, but may be 
more prevalent in structurally complex compounds such 
as YBCO. For example, YBCO's slightly orthorhom- 
bic structure facilitates crystal twinning, i.e. in a sin- 
gle crystal, there are generally domains with the nearly 
equal a and b directions interchanged. Separating such 
twin domains are well defined 45° grain boundaries or 
twin boundaries which have been shown to be effective 
extended vortex pinning sitesi^— Scanning tunneling 



microscopy (STM) imaging of a twinned YBCO crystal 
show that the areal vortex density is strongly modified 
by the twin boundaries?^ Small-angle neutron scattering 
(SANS) studies of YBCO confirm that the twin bound- 
aries strongly deform the VL. 15 ' 19 Understanding the in- 
fluence of such structural defects on the VL has been 
the subject of intense theoretical work,— ~— and is im- 
portant for two main reasons. Firstly, it affects p(B) and 
thus adds uncertainty to measurements of fundamental 
quantities like A and £, since it can be difficult to iso- 
late such extrinsic effects from changes in fundamental 
quantities of interest. Secondly, the degree of pinning of 
vortices determines the critical current density which is 
important for many applications. 23 

Measurements of the vortex state field distribution 
p(B) are most often done using SANS , 15 ' 19 nuclear mag- 
netic resonance (NMR),— and conventional muon-spin 
rotation (/iSR)^ All these methods probe the VL in the 
bulk and can be applied over a wide range of magnetic 
fields. It is also possible to probe the magnetic field dis- 
tribution p{B) near the surface of the sample using low 
energy-/iSR (LE-/iSR) in low magnetic fields. 2 - Recently 
we have demonstrated that similar information on p(B) 
near a surface can be obtained using /3-NMR.— This has 
the advantage that it can be applied over a wide range 
of magnetic fields. 

In this paper, we report measurements of the VL above 
the surface of the cuprate superconductor YBa2Cu307_a 
using /3-NMR. — ~— The 8 Li + beam was implanted into a 
thin silver overlayer evaporated onto several YBCO sam- 
ples. Measuring in the Ag allows one to isolate the con- 
tribution to p(B) from long wavelength disorder, i.e. dis- 
order that occurs on length scales much longer than the 
vortex spacing and A due to structural defects such as 
twin and grain boundaries. This is possible because the 
field distribution broadening just outside the supercon- 
ductor due to the VL inside has a very distinctive field 
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dependence. In particular, it vanishes in high magnetic 
fields where the VL spacing becomes less than the char- 
acteristic distance of the probe from the superconductor. 
On the other hand, long wavelength disorder has a much 
weaker dependence on magnetic field and dominates the 
observed p(B) in the high field limit. Our results show ev- 
idence for significant broadening of p(B) from such long 
wavelength disorder on the scale of D ss 1 fim, which is 
attributed to pinning at twin or other grain boundaries. 
The magnitude of the broadening is similar to that ob- 
served in bulk /iSR measurements, suggesting that the 
same broadening contributes to p(B) in bulk /iSR mea- 
surements. There is a crossover such that near Tc, where 
A > I), the broadening scales with the superfluid den- 
sity, whereas at lower temperatures, where A -C D, the 
broadening does not track the superfluid density. We 
discuss the consequences of this for the inference of A(T) 
from measurements of p{B) in polycrystalline supercon- 
ductors. 

The paper is organized as follows: section II reviews 
the theory for the field distribution, and its second mo- 
ment near the surface of a superconductor. Section III 
contains all the experimental details. In Section IV, we 
present the results. Finally in section V we discuss the 
results and draw conclusions. 



II. THE MAGNETIC FIELD DISTRIBUTION 
p{B) IN THE VORTEX STATE 



spacing given by the lattice constant of the crystal. Since 
this is much smaller than the VL constant, the resulting 
field distribution p(B) provides a random sampling of the 
spatially inhomogeneous field B(r) over the VL unit cell: 



p(B) = j J 5(B - B{r))dr 



(2) 



where the integral is over a unit cell of the VL of area 
A. In this paper we are concerned with the z-component 
of the magnetic field B z (parallel to the c-axis of YBCO 
samples), and refer to it simply as B. 

For an ideal triangular VL, the spatial dependence of 
the z-component of the magnetic field in or outside a type 
II superconductor follows the modified London equation, 



V 2 B 



(>2JJ B e(z) = |!©(z)5>-R). (3) 
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Here A = A a f, when the applied field is along the c-axis 
and the screening supercurrents flow in the ab plane, V 2 
is the 2D Laplacian, 8(z) = 1 for z > and zero oth- 
erwise, r is a 2D vector in the xy plane, R are the Bra- 
vais lattice vectors for the VL. We define the z axis as 
the normal to the surface of a superconducting slab with 
negative z outside the superconductor. The solution of 
Eq. is easily obtained using the Fourier transform, 
B(r, z) = Bo Ek e ' k z )i where the dimensionless 

Fourier components, -F(k, z), are given by^ 



In a type II superconductor, above the lower crit- 
ical field B c i, the magnetic field penetrates the sam- 
ple inhomogeneously forming a lattice of magnetic vor- 
tices, each carrying a flux quantum, $o = h/2e. In a 
perfect crystal, intervortex interactions lead to a long- 
range ordered 2D lattice of vortices, usually of triangular 
(hexadic) symmetry^ At the core of each vortex (a cylin- 
der of radius approximately the superconductor's coher- 
ence length £), the local magnetic field is maximal. Out- 
side the core, concentric circulating supercurrents par- 
tially screen the field which thus falls exponentially with 
a lengthscale A. The average magnetic field in the VL 
is the applied field Bq for flat samples, where demagne- 
tization effects are negligible i^2r— At a given field, the 
average vortex spacing, i.e. the lattice constant of the 
VL, a is fixed. For the triangular lattice this is 
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Here A 2 = k 2 + and k = ^[nx + ^^y] are the 
reciprocal lattice vectors of the triangular VL, where 
m, n = 0,±I,±2.... A cutoff function C(fc), approxi- 

mated by a simple Gaussian C(k) ~ e ~ , can be used 
to account for the finite size of the vortex core, where 
F(k,z) is replaced by F(k, z)C(k)£& However, the cor- 
rections due to C(k) are very small in our case, so it will 
be omitted. An approximate solution for the magnetic 
field along z (both inside and outside the superconduc- 
tor) is given by 



B(r, z) = B Y^ F(k, z) cos(k • r). 



(5) 
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If one considers the profile of the magnetic field along 
a line in the lattice (perpendicular to the direction of 
the applied field), it is thus corrugated with a period de- 
termined by a. This inhomogeneity in the magnetic field 
causes a characteristic broadening in local magnetic reso- 
nance probes such as the muons in /iSR or the host nuclei 
in NMR. Since the muon (or host nuclear spin) is at a 
well-defined lattice site(s), it samples the VL with a grid 



The second moment of p(B) at a depth z, a 2 — (B 2 ) — 
(B) 2 , where (..) is the spatial average, is given by 

a 2 =B 2 ]>> 2 (k,,). (6) 

The field distribution for a perfectly ordered triangu- 
lar VL calculated from Eqs. ([2]) and j5]), at B = 52 
mT and outside the superconductor at z = —90 nm and 
A = 150 nm (relevant to YBCO at T « T c ), is pre- 
sented in Fig. [TJa) (Ad = 0, defined below). It shows 
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FIG. 1: (Color online) (a) Simulation of p(B) in an applied 
field of 52 mT at a distance 90 nm from the superconductor 
using Eq. ([5]) convoluted with a Lorentzian of width Ad . (b) 
The broadening of p(B) versus the applied field. Solid lines 
represent a from Eq. © in the bulk, and 40 and 90 nm away 
from the superconductor. Long dashed line shows a at 90 nm 
from Eq. © for D — 4 (jm and / = 0.1. Inset: sketch of 
a possible vortex arrangement including vortex trapping at 
twin boundaries spaced by D and a regular triangular vortex 
lattice elsewhere, (c) cr(T) (normalized at T — 0), from Eq. 
© for B = 52 mT, z = -90 nm, and / = 0.1, is plotted 
against T/Tc for different values of D. A d-wave temperature 
dependence of A(T) is used.— In all figures A(0) = 150 nm. 



the characteristic high field skewness with a cutoff cor- 
responding to the field at the core of the vortices. The 
sharp peak corresponds to the most probable field B sa d 
at saddle points in B(r) midway between adjacent vor- 
tices. The low field cutoff occurs at the center of an 
elementary triangle of vortices. As we move farther from 
the superconductor, B sa d moves towards the applied field 
as the field approaches uniformity for z — > — oo, i.e., 
p{B) — > S(B — B ). This crossover occurs as exp(^z), 
as the z variation of the Fourier components FQc, z) in 
Eq. Q is controlled by k which takes values equal to or 
larger than 2-k /a. However, if instead we consider a sim- 
ple non-superconducting overlayer instead of free space, 
then the limiting p(B) will be the intrinsic lineshape in 
the overlayer material. 

As mentioned above p(B) is also affected by disorder 
in the VL due to pinning at structural defects in the 
crystal, where the superconducting order parameter is 
suppressed. Such disorder causes broadening of the mag- 
netic resonance, obscuring the features expected from an 
ideal VLj2& adding uncertainty to parameters of interest 
such as A and £. Relatively little is known about the de- 
tailed characteristics of this disorder. Accounting for the 
disorder of the VL is most often done by smearing the 
ideal lineshape with a Gaussian or Lorentzian distribu- 
tion of width Ad , where the latter is a phenomenological 
measure of the degree of disorderi^^ Calculated distri- 
butions for an applied field Bq = 52 mT are shown in 
Fig. nja) f° r = 0.3 and 0.6 mT, together with the 
ideally ordered VL (Ad = 0). Such disorder is more 
pronounced outside the superconductor and renders the 
lineshape symmetric when the depth dependent intrinsic 
VL broadening is smaller than Ad- 

One major difference between conventional /iSR and 
/3-NMR or LE-/1SR is the stopping range of the probe. 
In conventional /iSR the /i + stopping range is rs 150 
mg/cm 2 , yielding a fraction of a mm in YBCO. In con- 
trast, in /3-NMR or LE-/iSR, the mean depth of the probe 
can be controlled on a nm lengthscale from the surface. 
For implantation depths inside the superconductor, com- 
parable or larger than A, the /iSR lineshape (proportional 
to p{B)) is nearly field independent for 2B cl < B < B c2 
(for a < A), and the second moment of p(B) follows the 
formula^ 



0.00609$o 

A 2 (T) ' 



(7) 



neglecting the cutoff field. Using the latter makes a 
slightly field-dependent, but the corrections are small 
for fields Bq <C B C 2- Outside the superconductor, the 
magnetic field inhomogeneity of the VL vanishes over a 
lengthscale that depends on the spacing between vortices, 
a. In particular, the recovery to a uniform field occurs 
on a lengthscale of The field distribution is thus 

strongly field dependent when a (Bo) is of the order of 
\z\. This is shown in Fig. [ljb), where a due to the VL 
given in Eq. © is plotted against the applied field at a 
distance of 90 nm and 40 nm above the surface. In low 
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magnetic fields, the magnetic resonance lineshape out- 
side the superconductor is sensitive to both the intrinsic 
inhomogeneity of the VL as well as any additional broad- 
ening from disorder. However, in high magnetic fields the 
linewidth is dominated by VL disorder. 

Taking the view that the dominant source of disorder 
is due to twin or grain boundaries^ one can model the 
effect of disorder on the regular VL in different ways. 
The simplest is to assume that, in addition to the regular 
triangular lattice, a fraction of vortices is trapped along 
the structural defects such as twin or grain boundaries 
as shown in the inset of Fig. [lib). The local field in real 
space will be the superposition of both contributions 



B{v,z) = B v[ (r,z) + B dls (r,z), 

= < ]T F(k, z)e 4k r + < s £ F(g, zje*', 



(8) 



where B vl (r 7 z) is the field due to the regular VL, 
B dis (r, z) the field due to the vortices pinned by disorder 
and g is some generally incommensurate wave vector re- 
lated to the pinning, which for simplicity we take to be 
of the form g = ^nx + ^roy, where a' is the spacing 
between vortices in a boundary, and D is the separa- 
tion between boundaries as drawn in the inset of Fig. 
IHb). The average field is then B = VW) 2 + {B dis ) 2 , 
where B dls = fBo and / is the fraction of pinned vortices 
(0 < / < 1). The second moment of B(r, z) from Eq. flS)) 
can be then easily calculated: 



scales with 1/A 2 as predicted for an ideal VL (see Eq. 
([7])). In particular, at low T, the broadening is almost 
T-independent irrespective of superconducting gap struc- 
ture. It is interesting to note that the first /iSK studies 
on powder samples of cuprates showed a very flat vari- 
ation in the linewidth j 41 ' 42 This was taken as evidence 
for s-wave superconductivity. Later measurements on 
high quality crystals of YL^CusOy-a showed a much 
different low temperature behaviour^ and, in particu- 
lar, a linear variation in 1/A 2 (T) consistent with d-wave 
pairing. 44 Although the lineshapes in powders are ex- 
pected to be more symmetric than in crystals due to the 
additional disorder and random orientation, the differ- 
ent temperature dependence is surprising since it was 
thought that the line broadening from disorder should 
also scale with 1/A 2 (T)£5. The current work provides a 
clear explanation for the discrepancy between powders 
and crystals. In powders, the line broadening is dom- 
inated by long wavelength pinning of vortices at grain 
boundaries. Consequently the resulting broadening at 
low temperature reflects variations in the vortex density 
and is thus only weakly dependent on temperature. In 
later work on crystals, the contribution from such long 
wavelength pinning is much less important. This is ev- 
ident from bulk juSR in crystals where one observes the 
expected characteristic lineshape associated with a VLi2£ 



III. EXPERIMENTAL DETAILS 



(l-/ 2 )£F 2 (k,z) + / 2 ][> 2 (g,z) 

k#0 g^O 



(9) 



It is clear from Eqs. (Ql and ©, that the broadening 
from the VL (first term) at a distance z outside the super- 
conductor becomes small at high magnetic fields where 
\z\ a. However the broadening outside the supercon- 
ductor due to disorder (second term) remains large pro- 
vided | z | is not much larger than Since D and / de- 
pend on the arrangement of twin boundaries we expect 
them to be sample dependent. In addition, one may also 
anticipate that / will decrease at high magnetic fields 
where the increased repulsive interaction between vor- 
tices overcomes vortex-pinning. Therefore, we assume a 
simplified phenomenological parameterization / = 5B Q 7 , 
where S is temperature and sample dependent and 7 > 1. 

The broadening of the field distribution due to a reg- 
ular VL can be significantly larger when introducing the 
effect of disorder due to the twin and grain boundaries. 
When taking the disorder into account, a of Eq. © is no 
longer zero at high magnetic fields as seen in Fig. [TJb). 
This is because the broadening has a disorder compo- 
nent which decays on a length scale of D rather than 
a, where D ^> a (we also assume D ^> a' and thus ig- 
nore the effect of the spacing within the twin bound- 
aries). Consequently a shows a strong deviation from 
the ideal VL result as seen in Fig. [T}c as D increases. 
In this case the second moment from Eq. @ no longer 



The measurements were carried out on three differ- 
ent near-optimally doped YBCO samples, two flux-grown 
single crystals and a thin film. I) The twinned single crys- 
tal in the form of a platelet ~ 0.5 mm thick with an area 
~2x3 mm 2 had Tc = 92.5 K. It was mechanically pol- 
ished with 0.05 /im alumina, then chemically etched with 
a dilute (0.8%) Bromine solution followed by annealing 
at 200° C in dry N2 to improve the surface quality. It 
was then sputter coated with a 120 nm thick Ag film 
(99.99 % purity) at room temperature in an Ar pressure 
of 30 mTorr. The deposition rate was 0.5 A/s, and to 
ensure Ag uniformity, the crystal was rotated. II) The 
optimally doped detwinned single crystal had Tc = 92.5 
K, ~ 0.5 mm thickness, and area ^3x3 mm 2 . The crys- 
tal was cleaned, annealed, and mechanically detwinned. 
A 120 nm thick Ag, from the same source as above, was 
sputtered onto the prepared surface under similar con- 
ditions. Ill) The film of Tc = 87.5 K, critical current 
density Jc — 2.10 6 A/cm 3 and 600 nm thickness, sup- 
plied by THEVA (Ismaning, Germany), was grown by 
thermal co-evaporation on a LaA10*3 substrate of area 
9x8 mm 2 . The film was coated in situ with a 60 nm 
silver layer (99.99% purity). 

The experiments were performed using the /3-NMR 
spectrometer at the ISAC facility in TRIUMF, Canada, 
where a highly nuclear-spin-polarized beam (intensity 
~ 10 6 ions/s) of 8 Li + is produced using collinear opti- 
cal pumping with circularly polarized laser lights The 
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FIG. 2: (Color online) Implantation profiles of Li at ener- 
gies of (a) 5 keV into 60 nm of Ag with the mean at 40 nm 
away from YBCO film, and (b) 8 keV into 120 nm of Ag with 
the mean at 90 nm away from YBCO crystals, as calculated 
via TRIM.SP. 45 Solid lines are phenomenological fits. 



beam is directed onto the sample which is mounted on the 
cold finger of a He flow cryostat and positioned in the cen- 
tre of a high homogeneity superconducting solenoid. The 
beamline and entire spectrometer are maintained in ul- 
trahigh vacuum (ICP 10 Torr). In /3-NMR measurements, 
the 8 Li + nuclear spin polarization is monitored via its 
asymmetric radioactive beta decay (lifetime r = 1.203 
s), where the high energy (several MeV) beta electron is 
emitted preferentially opposite to the nuclear spin direc- 
tion. The experimental asymmetry, defined as the ratio 
p^_g of the count rates in two plastic scintillation detec- 
tors placed in front (F) and at the back (B) of the sample, 
is proportional to the probe's spin polarization . 27 ' 28 

The whole spectrometer can be biased at high voltage, 
allowing one to tune the implantation energy of 8 Li + ions 
and their implantation depth between 5-200 nm. There- 
fore, the implanted 8 Li + can monitor the depth depen- 
dence of the local magnetic field distribution in materials 
at nm scale by measuring the NMR lineshape in a man- 
ner analogous to conventional NMR^ i 28 ' 29 In this work, 
the 8 Li + ions are decelerated to stop in the Ag overlayer 
deposited on each of the three YBCO samples. Implanta- 
tion profiles of 8 Li + were calculated using the TRIM.SP 
code^ examples of which are shown in Fig. [3] The 
implantation energies used in this study (8 keV in the 
crystal samples and 5 keV the film), were tuned to stop 
all the 8 Li + within the Ag. The mean distances are 90 
and 40 nm from the Ag/YBCO interface in the crystals 
and film, respectively. 

The /3-NMR measurement is carried out by monitor- 
ing the time averaged nuclear polarization through the 



beta decay asymmetry, as a function of the radio fre- 
quency (RF) co of a small transverse oscillating magnetic 
field Bi = B 1 cos(ut)x, where B 1 ~ 0.01 mT. When w 
matches the Larmor frequency WLi = TLi-Blocalj where 
7Li = 6.3015 kHz/mT is the gyromagnetic ratio and 
-Biocai is the local field, the 8 Li + spins precess about 
Siocai, causing a loss of polarization. To establish the 
vortex state in the YBCO samples, they are cooled in a 
static magnetic fields Bo > B c i applied parallel to the o 
axis of YBCO (normal to the film and platelet crystals). 
Bq is also parallel to both the initial nuclear spin polar- 
ization and the beam direction. The local field sensed by 
the 8 Li + is determined by the applied field and the in- 
ternal magnetic field generated by the screening currents 
associated with the vortex lattice. Thus, -Bi OC ai is dis- 
tributed over a range of values, which can be calculated 
using 

P (B) = J dzp(z)jj &r8(B-B(r,z)). (10) 

where p{z) is the implantation profile calculated using 
TRIM.SP given in Fig. H 

When 8 Li + is implanted in Ag (with no superconduct- 
ing substrate) at temperatures below 100 K, it exhibits 
a single narrow resonance at the Larmor frequency^ 
The resonance should yield an approximately Gaussian 
distribution caused by nuclear dipolar moments. 46 How- 
ever, continuous wave RF leads to a power-broadened 
Lorentzian lineshape, whose linewidth is small (1 kHz 
w 0.15 mT) and corresponds to the dipolar broaden- 
ing due to the 107 > 109 Ag nuclear moments and RF power 
broadening.— In the presence of any additional magnetic 
inhomogeneity in the Ag, due for example to a VL as- 
sociated with a superconducting substrate, the observed 
resonance lineshape will be a convolution of the narrow 
RF power broadened Lorentzian of Ag with the (depth 
dependent) field distribution due to the VL in the sub- 
strate. There are unique aspects of measuring the field 
distribution in the Ag overlayer compared to the super- 
conductor itself. As mentioned above, in high magnetic 
field, it is possible to isolate and study the broadening 
due to VL disorder that occurs on a long length scale. 
Also in low magnetic fields, where the broadening is dom- 
inated by the VL, it should be possible to measure A 
in magnetic superconductors since the field distribution 
above the sample is free of any internal hyperfine fields 
that make a bulk measurement impossible.— 



IV. RESULTS 

The /3-NMR resonances were measured as a function of 
temperature under field-cooled conditions at fields rang- 
ing from Bo = 20 mT to 3.3 T in each one of the three 
samples. Fig. [31 shows typical resonance lineshapes at 
various temperatures in sample I with Bo = 51.7 mT. 
Above Tc, the line broadening is small and temperature 
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FIG. 3: (Color online) /3NMR resonances in Ag/YBCO (crys- 
tal I) at temperatures 100 K, 80 K, 20 K, and 4.5 K measured 
in a magnetic field Bo of 52.3 mT applied along YBCO c-axis. 
The solid lines are best fits using a Lorentzian. 



independent as expected from nuclear dipolar broaden- 
ing. Below Tq, the field distribution in the Ag overlayer 
broadens dramatically from the VL in the underlying su- 
perconductor. Such broadening was observed in all sam- 
ples and at all magnetic fields, although there are sig- 
nificant variations as a function of both magnetic field 
and sample as discussed below. The first thing to note 
is that the lineshape is very symmetric and fits well to a 
simple Lorentzian. This is much different from the asym- 
metric lineshape observed with conventional /j,SR in sam- 
ples similar to I and IL 25 i 43 The other significant differ- 
ence between the current results and previous bulk fiSK 
measurement o 4 ' 9 ' 25 i 43 on crystals is that the broadening 
at low temperatures is only weakly dependent on temper- 
ature, as may be seen by comparing the resonances at 20 
K and 4.5 K. In contrast, the broadening from an ordered 
VL lattice scales with 1/A 2 and consequently in YBCO 
shows a strong linear T-dependence at low temperatures 
due to the ri-wave superconducting orderi 43 ' 44 

The observed lineshape fits well to a convolution of 
two Lorentzians, one from vortices in the superconduct- 
ing state with a full width at half maximum (FWHM) 
A sc , and one from other sources determined from the 
normal state of FWHM A ns . The width of a convolution 
of two Lorentzians is the sum of the individual widths: 
A(T) = A SC (T) + A ns . Therefore, the contribution from 
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FIG. 4: (Color online) The vortex-related broadening below 
T c , A SC (T) = A(T) - A ns of the twinned (full symbols) and 
detwinned (opaque squares) YBCO crystals in an applied field 
Bo- A(T) is the linewidth at temperature T of the Lorentzian 
fits and A ns is the constant linewidth in the normal state. 
Solid lines represent a fit using Advl = 2.355<r where a is 
given in Eq. ([9]) and D and / are varied to fit the data. A 
d-wave temperature dependence of X(T) in YBCO is used, 37 
where A(0) = 150 nm. 



vortices in the superconducting state can be obtained by 
simply subtracting the temperature independent normal 
state width. Fig. Q] shows the resulting A SC (T) as one 
enters the superconducting state in samples I and II. At 
low field, the measured width (~ 2.2 mT) at low tem- 
perature is larger than expected from a regular VL and 
decreases significantly in the detwinned crystal to about 
~ 0.6 mT. For comparison, simulations using Eq. (|10[) 
and the 8 Li + stopping profile in Fig. &,b); indicate that 
the broadening due to a regular VL is only Avl ~ 0.3 
mT. At 3.33 T, the discrepancy between the observed 
width (see Fig. 0} and that expected from a regular VL 
is even more dramatic. At this high field the vortices are 
spaced so closely (a s» 27 nm), that there should be no 
detectable broadening from a regular VL for our stopping 
depths. This can be seen clearly from the simulation in 
Fig. [ljb), where the VL broadening approaches zero at 
high fields. In contrast, the data at 3.33 T shows sig- 
nificant broadening below Tq which is therefore solely 
attributed to vortex disorder on a long length scale. 

The temperature dependence of the broadening is also 
much weaker than expected from a regular VL in YBCO, 
where 1/A 2 has a strong linear term due to the d-wave 
order parameterj 43 ' 44 The observed temperature depen- 
dence fits well to our model of disorder, where A SC (T) 
is compared to an estimate of the FWHM given by 
Advl ~ 2.355er, where a is given in Eq. ©. This leads 
to an estimate of D of the order of a micron, consistent 
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FIG. 5: (Color online) Comparison of the field distributions in 
the three samples taken at temperatures 100 K (top panel), 5 
K (crystals) and f0 K (film). The x-axis is shifted by Bo, the 
applied field which is 52.3 (51.7) mT for the crystals (film). 
Solid lines are Lorentzian fits and dashed lines are the simu- 
lation described in the text. 
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FIG. 6: (Color online) Superconducting broadening A SC (T) 
of the /3-NMR resonance spectra at temperatures ~ 4.5-10 K 
in the three samples. The experimental broadening is com- 
pared with the linewidth of an ideal VL Avl (dashed lines), 
and Advl of a disordered VL (solid lines), which are both 
weighted by the 8 Li + profile given in Fig. [2] 



with the the separation between twin boundaries or grain 
boundaries^ 8 - In the detwinned crystal, D is found to be 
larger but not infinite since the detwinning is not com- 
plete. The fraction of vortices / pinned by structural 
defects in the twinned crystal is about ^0.1 at low field 
(52 mT) and decreases considerably at high field (3.33 
T). Thus, the amplitude of the enhanced vortex density 
at the defects, fBo, varies between 2-4 mT at all fields. 
In the detwinned crystal, / ~ 0.01, is an order of mag- 
nitude smaller than in the twinned crystal at the same 
field, with small variation in the vortex density (0.5 mT) 
compared to the twinned crystal. These results are con- 
sistent with expectations from pinning at twin bound- 
aries. In particular, one expects the fraction of vortices 
pinned will decrease in the partially detwinned crystals. 
Also, it is reasonable to expect that in a high magnetic 
fields the fraction of vortices pinned will decrease due 
to the smaller separation between vortices and resulting 
increase in the repulsive interaction. 

In Fig. [SJ the spectra in all three samples above and 
below Tq are compared with the corresponding simulated 
field distributions. The observed lineshapes are all sym- 
metric, and significantly broader than expected, showing 
little or no sign of the characteristic VL field distribu- 
tion. Simulation of the VL lineshapes (dashed lines) was 
done using Eqs. (gj, (J5J), and (fTU)) . for A = 150 nm, and 



was convoluted with a Lorentzian representing the nor- 
mal state spectra with A ns = 0.3 mT. The theoretical 
lineshape for the film is broader and asymmetric because 
it is weighted by the 8 Li + stopping distribution which was 
on average closer to the superconductor. The lineshapes 
for the crystals are almost symmetric as the ideal VL line- 
shape at the depths of an average 90 nm away from the 
superconductor are narrower than the Lorentzian they 
are convoluted with. The magnetic field dependence of 
the superconducting linewidth A SC (T) at low tempera- 
tures ~ 4.5-10 K is plotted in Fig. [5] In all samples, we 
find that the broadening is largest at low field and de- 
creases gradually with increasing field. Also, in all cases 
the broadening remains large and well above the predic- 
tion from a regular VL, approximated by Avl ~ 2.355<r 
where a of an ideal VL is given in Eq. ([6]) and weighted by 
the 8 Li + profile given in Fig. [5] The broadening is sub- 
stantially smaller in the detwinned crystal compared to 
the other samples. One can account for all of the data us- 
ing a linewidth due to a disordered VL, Advl ~ 2.355c, 
where a is now given in Eq. ^) and weighted by the 
8 Li + profile plotted in Fig. O The data is well fitted 
(see Fig. [5]) by assuming that the twin/grain boundaries 
spacing D is sample dependent of the order of a few mi- 
crons, and by assuming a phenomenological form for the 
fraction of pinned vortices / = SBq' 1 with 7 ~ 1.1 and S 
sample dependent. 
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V. DISCUSSION AND CONCLUSIONS 

It is clear that the /3-NMR lineshapes in the Ag over- 
layer differ substantially from that expected from a well- 
ordered VL field distribution. This has little to do with 
the method of observation. For example, in the con- 
ventional superconductor NbSe 2 , /3-NMR shows the ex- 
pected VL lineshape.— The lineshapes reported here in 
the YBCO film are also qualitatively different than that 
seen with LE-/1SR in a YBCO film coated with a 60 nm 
thick Ag layer.— In that experiment the authors found 
a more asymmetric lineshape in the Ag overlayer which 
was closer to that of a regular VL. Some of this differ- 
ence may be due to the different pinning characteristics 
of the samples, although the YBCO film used by Nie- 
dermayer et al. was from the same source as sample III. 
Also, the LE-/iSR experiment was probing the VL closer 
to the interface and in a lower applied field where the 
disorder is less important compared to the contribution 
from the ordered VL. The symmetry and large broaden- 
ing of the lineshape at low fields cannot be accounted for 
by VL melting (at a reentrant vortex liquid state near 
B c \) which would instead yield a motional narrowing of 
the field distribution^ 

The observed resonances in the current experiment are 
dominated by long range variations of the vortex density 
across the face of the sampled Such disorder in the VL 
can produce a symmetric lineshape ; 38 ' 50 and can broaden 
the field distribution significantly compared to that of the 
corresponding ordered stated Weak random pinning or 
point-like disorder due to oxygen deficiency may slightly 
distort the VL, and may also broaden the lineshape^ 
However, the correlated disorder due to the twin and 
grain boundaries is dominant at long wavelengths ; 40 ' 52 
and therefore we are mostly sensitive to the twin/grain 
boundaries. Indeed, the position of the probe outside 
the superconductor enhances its sensitivity to long wave- 
length disorder, as the proximal fields fall off with dis- 
tance as exp(— 2ir\z\/D) where D is the wavelength of 
the inhomogeneity of the field M- The broadening is re- 
duced in a detwinned crystal where the twin boundaries 
are more sparse as shown in Fig. [SJ thus the vortex den- 
sity variation across the face of the sample is smaller than 
in the twinned crystals as / is largely reduced. 

The extrinsic broadening due to disorder at low tem- 
perature, Ad = A sc — Avl, reported here is between 
0.5 mT and 2.5 mT. This is remarkably close to the ad- 
ditional Gaussian broadening required to explain line- 
shapes in bulk /iSR measurements on crystals! 25 ' 38 ' 54 In 
the bulk, this extrinsic broadening is small compared to 
the intrinsic VL broadening, whereas outside the sam- 
ple the reverse is true. It is important to note the 
T-dependence of the extrinsic broadening in Fig. @] 
does not follow the superfluid density (cx 1/A 2 ) which 



varies linearly at low-T because of the ci-wave order 
parameter! 44 ' 55 Instead, we observe a much weaker In- 
dependence. This is expected from our model of disor- 
der which occurs on a long length scale D. For exam- 
ple, at low-T where A is short compared to D, the flux 
density outside the sample is determined solely by in- 
homogeneities in the vortex density, and is independent 
of the superfluid density as the vortices are static and 
well-pinned in the twin boundaries. The current results 
may also explain early ^iSR work on HTSC powders and 
sintered samples which mistakenly indicated an s-wave 
T-dependence of 1/A 2 i 41 ' 42 It is likely in these cases the 
linewidth was dominated by extrinsic VL disorder on a 
long length scale. This tends to flatten the T-dependence 
of the linewidth and the effective A obtained from the 
analysis ! 41 ' 56 Therefore we conclude that although the 
linewidth obtained from powders can be useful in mak- 
ing rough estimates of A, one cannot extract accurate 
measurements of A or its T-dependence without addi- 
tional information about the source of broadening and in 
particular VL disorder. 

In conclusion, we have measured the magnetic field 
distributions due to the vortex state of YBCO using (3- 
NMR. We find a significant inhomogeneous broadening 
of the NMR attributed to the underlying VL in YBCO. 
However, the observed resonances have several unex- 
pected properties. In particular, they are broader and 
more symmetric than for an ideal VL. The anomalous 
broadening is most evident in high field where there is 
no significant contribution from the regular VL. These 
effects are attributed to long wavelength disorder from 
pinning at twin or grain boundaries. The temperature 
dependence of the disorder-related broadening does not 
scale with 1/A 2 , suggesting there is a contribution to the 
linewidth in the bulk of the vortex state that does not 
track the superfluid density. This is likely to have only 
a minor effect on the interpretation of data on crystals 
where the observed lineshape is close to that expected 
from a well ordered VL. However, it can be significant in 
powders or crystals where there is substantial disorder in 
the VL. In particular, when the broadening is dominated 
by VL disorder on a long length scale (i.e. much big- 
ger than A) the temperature dependence of the linewidth 
does not scale with 1/A 2 and therefore cannot be used to 
determine the symmetry of the superconducting gap. 
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